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ABSTRACT 

Experiment E865 @) at the Brookhaven AGS was set up primarily to search for 
the lepton flavour violating decay K + — ► 7r + /i + e~ {K niie ) with high sensitivity. 
The flexibility of the apparatus allowed also to obtain more than an order 
of magnitude larger than previously available event samples on the following 
decay modes: 7r+e+e~ (K^ ee ), tt + (j, + ij,~ {K niltl ), -K + -K~e + v e (K ei ), ^+e+e~z/ p , 
and e + e + e~v e . The report focusses on the K nfJ , e results and those on other 
lepton flavour violating decays as well as the K e ± data, from which a new, quite 
precise value for the s-wave irir scattering length can be deduced. 



1 The E865 experiment: goals and techniques 

The search for a violation of the apparent, but theoretically not well founded 
lepton number conservation law continues at several laboratories around the 



world. Improving the sensitivity beyond existing limits for the decays [i + — > 
e+7 (< 1.2 x 10" 11 ) I) and /x" A -> e - A' (< 1.7 x 10~ 12 ) I) is the goal of ongo- 
ing or planned experiments at the Paul Scherrer Institut (PSI) and Brookhaven 
National Laboratory (BNL) @) . The observed evidence for neutrino oscilla- 
tions H) has further stimulated the interest in this area. Since extensions of 
the Standard Model like supersymmetry, leptoquarks, technicolor and other 
horizontal gauge boson models have been proposed which allow the possibility 
of flavour changing neutral currents and lepton family violating interactions, 
the forbidden K decays K° L -> A* ± e =F (< 4.7 x lCT 12 )*!) and K+ -> TT+fi+e" 
[K-niie < 2.0 x 10~ 10 ) llP were also investigated. E865 was set up to search for 
the latter decay with improved sensitivity. 




Figure 1: Plan view of the E865 spectrometer. 



The E865 detector (see Fig. P resided in an unseparated 6 GeV/c beam 
of 1.5 x 10 7 K + accompanied by about 3 x 10 s ir + and p per AGS cycle. About 
8% of the K + decayed in the 5 m long decay volume. The decay products were 
separated by charge by a first dipole magnet and then momentum analysed by a 
second dipole magnet sandwiched between two pairs of proportional chambers. 
The wire chambers with four planes each were kept at lower voltage in the region 
were the beam passes. This arrangement yields a resolution of cr p /p 2 ~ 0.003, 
where the momentum p of the decay products (7r ± ,/i ± or e ± ) ranges between 
0.6 and 3.5 GeV/c. Four gas filled Gerenkov counters were positioned inside 
and after the spectrometer magnet for e identification. An electromagnetic 



calorimeter of Shashlyk type u) further helped to separate e from other decay 
products. It consists of 600 modules, each 15 radiation lengths deep in beam 
direction. It is followed by an array of 12 planes of proportional tubes separated 
by iron plates used to discriminate pions against muons. 

The event trigger is based on the information from four scintillator ho- 
doscopes, D and A before and after the spectrometer, B and C in front and in 
the middle of the muon stack. In the lowest trigger level three charged particles 
are required, each identified by a cluster in the calorimeter and a signal from 
the A-counter in front of it. Depending on the final state the particle identifi- 
cation signals and topological information is added at higher levels, e.g. a K^^e 
candidate is given by an e~ (Cerenkov on beam left filled with H2 at atmo- 
spheric pressure) and two particles on the beam right, one traversing the muon 
stack, vetoed with signals from the CH4 filled Cerenkov on beam right. The 
flexibility of the trigger and particle identification system allowed to accumu- 
late high statistics data samples either concurrently with K nfie high intensity 
running or in dedicated runs at lower intensity on the decays 7r + e + e~ (K vee ), 
7T + n + fj,~ (Kirfift), ■n + ir~e + v e {K e ^), /i + e + e~f M , e + e + e~v e for physics analysis 
and on 7r + 7r° (e + e~7) (KDaiitz) and 7r + 7r + 7r _ (K T ) for normalisation. 

The report focusses on the K n/ie results and those on other lepton flavour 
violating decays as well as the K e ± data, from which a new, quite precise value 
for the s-wave nir scattering length can be deduced. This quantity plays a 
central role in the low energv effective theory of strong interactions - chiral 
perturbation theory (ChPT) Both the flavour changing neutral cur- 

rent decays K vjl(L and K^ ee , which will not be discussed here in detail, are 
known to be dominated in the Standard Model by long-distance effects in- 
volving one-photon exchange, and are of interest for ChPT tests, too. Our 
published results 0' for these decays include precise branching ratios and 
form factor information, establish firmly the vector nature of the interaction, 
resolve the discrepancy of older data with efj, universality, and provide evidence 
for the existence of a pion loop term calculated in next-to-leading order ChPT. 

2 Lepton flavour number violating decays 

Not all our K vlLe data have been analysed yet. The data taken during the last 
running period in 1998 are projected to push the sensitivity limit below the 
10~ n branching ratio level. Table |l| lists the limits deduced from the 1995 El) 



Table 1: Upper limits for lepton flavour number violating decays established 
by experiment E865 : LFNV: violation of lepton flavour number con- 

servation (L^, L e ); LNV: violation of total lepton number conservation (L = 

L^ + L e ); Gi M . number of quarks and leptons of generation i ; GV: violation 
of generation number conservation. (PDG: particle data group IZ^ ) 
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Branching 
ratio 
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and 1996 data 0) . Figure || shows the mass distribution of the few K niJ , e 
candidates remaining after a likelihood analysis. No event is found with a like- 
lihood for the K^^e hypothesis exceeding 20%. Background events may come 
from K + decays with errors in kinematic reconstruction or particle identifica- 
tion measurements or from accidental combinations of three tracks (dominant 
contribution). The probabilities to misidentify a n~ as an e~ , or an a 7r + as an 
e+ are 2.6 x 10~ 6 and 1.7 x 10~ 5 (0.049), while the identification efficien- 
cies for 7r + , ix + and e~ are 78%, 74% and 55%, respectively. The likelihood 
distributions constructed for these background events let us then expect 2.6 
events below the 20% line indicated in Fig. || consistent with the threeevents 
observed. Combined with earlier data @) our results yield an upper limit for 



the Ki^e branching ratio of 2.8 x 10~ n (sec Table Q). 

In models described by a horizontal gauge interaction EH' EH) one may 
use this branching ratio to deduce a limit for the mass of an intermediate boson 
X by comparison to the allowed decay K + — ► 7r°^ + ^ Al (K^): 



B.R.(K^e) 16 fgxVfMwY 



B.R.{K^) ~ sin 2 f? c \g ) \M X J 



{g/gx)M x < 60 TeV/c 2 



Here 9 C is the Cabbibo angle, and gx and g are the couplings for the new and 
the weak interaction, respectively. 

Unlike K V n ei which only violates lepton flavour conservation, the other 
decays for which Table [I] lists upper limits, also violate generation number 
conservation ED . We used our K^^^ and K e ± data samples to search for these 
decays. These data sets were taken with dedicated triggers at lower beam 
intensity. The and K efi7! mass plots are shown in Fig. [2] as examples. 

The most significant reduction of the existing limit El) (by a factor of 50000) 
was obtained for K^^, a decay which is sensitive to the same exotic mechanism, 
which may permit nuclear double /3-decay 113 . 

3 K e 4 and the tttt scattering length 

Since more than 30 years, the K e ± decay has been used to study 7T7r-scattering 
at low energies. The most recent experiment performed by a Geneva-Saclay 



collaboration at CERN accumulated 30000 events. Using these data and 
some peripheral ttN — > irirN data did lead to the presently generally accepted 
value for the isoscalar s-wave scattering length of a° = (0.26 ± 0.05) (in units 
ofM-i). 

On the theoretical side tttt scattering is the golden reaction for chiral per- 
turbation theory (ChPT) © . This low -energy effective theory for the strong 
interaction makes firm predictions for the scattering length. The pioneering 
tree level calculation by Weinberg E3) yielded a° = 0.156 (a 2 , = —0.046). The 
one-loop (a[j = 0.201, a 2 , = —0.042 0) and the recently completed two loop 
calculation (a[j = 0.217, a 2 , = —0.041 IH' EIP) show a satisfactory convergence. 
The most recent calculation matches the known chiral perturbation theory 
representation of the tttt scattering amplitude to two loops with a phenomeno- 
logical description that relies on the Roy equations (see below) . The corrections 
to Weinberg's low energy theorems for the s-wave scattering lengths are worked 




Figure 2: Top: scatter plot of log-likelihood versus invariant mass for K vjle 
candidate and Monte Carlo generated signal events; the horizontal lines demark 
the 3(7 mass region 0). Bottom right: as above with coordinates exchanged 
for K efJ-7r candidate events; the box indicates the signal region. Bottom left: 
invariant mass distribution for K,.,.^ candidate events; the dashed line indicates 

the expected signal \L31 . 



out to second order in the expansion in powers of the quark masses, and the 
uncertainties from higher order effects are estimated reliably for the first time. 
The resulting predictions = 0.220 ± 0.005, and ag = -0.0444 ± 0.0010 are in 
remarkable agreement with our new experimental result, which we will discuss 
in the following. 

Figure shows these predictions and also a number of older, alternative 
model theoretical predictions . Since the extraction of the irir amplitude in 
the high statistics irN — > irirN experiments is model dependent, and the irir 
atom data from the DIRAC experiment are not yet available, one had to 
rely primarily on the K e 4 data for the comparison between theory and experi- 
ment, which were not sufficiently accurate to discriminate between the different 
models. 



[x 100] 




Figure 3: Left: previous measurements and predictions for the scattering 
length (taken from Ref. El) ). Right: azimuthal asymmetry averaged over all 
other kinematical quantities with standard phase irir shifts and with phase shifts 
set to zero. Insert: kinematical quantities used in the description of decay. 



K e 4 events are selected with a vertex within the decay tank of acceptable 
quality, a momentum reconstructed from the three daughter particles below 
the beam momentum. The momentum vector of the incoming K + is obtained 
from the hit in the pixel counter (see Fig. and the vertex, and requiring 
the K + to come from the production target 27.5 m upstream as well as an 
unambiguous identification of the e + and the ir~ in both Cerenkov counters 



and the calorimeter. The event sample contains a small amount of background 
events, mainly from K T decay (1.3 ± 0.3 %) and accidentals (2.6 ± 0.2 %). A 
K T can fake a K e ^ by either a misidentification of one of the 7r + due to J-rays, 
photomuliplier noise or the presence of an additional positron, or a decay of 
one 7r + directly or via a /i + into a e + . The dominating accidental background 
is a combination of a n + Tt~ pair from K T decay with a positron from either 
the beam or from a K^ai decay. 

After event selection 406T03 events remain, of which we estimate 388270 
± 5025 to be K e ± events. This corresponds to an increase in statistics by more 
than a factor 12 compared with previous experiments E^). 

Quantitative understanding and analysis of the data relies strongly on 
the quality of Monte Carlo simulation of the detector response. This GEANT- 
based HH) simulation is supplemented with separately measured efficiencies of 
the various detector parts, and was extensively tested against data in the vari- 
ous decay channel under study. Figure ^ shows a few control plots demonstrat- 
ing the good agreement between measured and simulated distributions of the 
kinematic variables. The total Monte Carlo data sample comprises 81.6T0 6 
generated K e 4 events of which 2. 9 TO 6 were accepted after reconstruction and 
passage through the same analysis codes as the data. We also compare Monte 
Carlo with data distributions for the kinematically very distinct K T and Kdai 
decays, and are able to reproduce the measured branching ratios for these de- 
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cays. A chiral perturbation theory calculation on the one loop level 
is used to model the physics of the decay at this point. Radiative corrections 
are included following Diamant 

For the description of the K e 4 decay three reference frames are commonly 
used (see Fig. ||): the K + rest system (£k), the n + n~ rest system (£„-„•) and 
the e + v rest system (E ev ). The decay kinematic is then fully described by five 
variables: s T = M£ w , the invariant mass squared of the dipion; s e = M^ v , the 
invariant mass squared of the dilepton; 9 n , the angle of the n + in with 
respect to the direction of flight of the dipion in E^; 9 e , the angle of the e + in 
£ e „ with respect to the direction of flight of the dilepton in <j>, the angle 
between the plane formed by the pions in E W7r and the corresponding plane 
formed by the leptons. The experimental resolution for these five variables is 
cr(s„) = 0.00133 GcV 2 , o-(s e ) = 0.00361 GeV 2 , a(6» 7r ) = 0.147, a(9 e ) = 0.111 
and a((f>) = 0.404. 



The K e 4 branching ratio is determined relative to K T decay (Br(r) = 
5.59 ± 0.05 %), for which events were collected with a minmum bias trigger 
concurrently with K e 4 events. We obtain 



BR^) = (4.109 ± 0.008 (stat.) ± 0.110 (syst.)) • 10~ 5 

The result is in agreement with the average of previous measurements: (3.91 ± 
0.17) • 10 -5 ll^. The largest of the 17 contributions to the systematic er- 
ror (2.68%), which we have evaluated, stem from the background subtraction 
(1.2%), the Cerenkov efficiencies (1.5%) and the K T branching ratio (0.9%). 

The analysis of the K e ^ decay distributions is far from trivial, and requires 
a fair amount of theoretical input. We list below the basic ingredients and 
prepare the notation for the subsequent result section. 

The most general form of the decay matrixelement is written most conve- 
niently using the following combinations of four-momenta and Lorentz invari- 
ants: 

P=Pl+P2, Q=Pl-P2, L = <p e + Pv , X= [(P ■ L) 2 - SjrSe] 1 / 2 , 

M - ^V: s u(pM1 - J 5 )v(Pe)(V» - A^) , (!) 

Afl = -M- ( Fpt * + G ^ + RL>1 ) > yiL = -£r efl " paL » p pQ° ■ 

F, G, R, H dimensionless complex functions of Px • P2, Px ' P, P2 ■ P, or of s v , 
s e and 9 n . If terms cx M 2 /s e are neglected R does not contribute to the decay 
rate, first evaluated by Pais and Treiman HI) : 
Q2 y2 2M n q 

dT 5 = 2 i2^6^5 ^ 7= J5(s-K,s e ,e„,6 e ,(f))ds v ds e dcos9„dco86 e d^ . (2) 



For the detailed decomposition of the factor J5 in the above expression in terms 



of the kinematical variables we must refer to Ref. — 1 E^) . A convenient partial 
wave expansion of the form factors in the variable 9^ gives the following 0) : 

+f P {q 2 /s,) 1/2 (P ■ L) cos^ exp(^K^)) , 

(3) 

G = {g P + g' p q 2 + g e s e ) exp(i5{(s„)) , 
H = (h p + h' p q 2 )exp(i8l(s n )) , 

where q is the pion momentum in E TTr , divided by M^. This parameterization 
yields ten new form factors f s , f' s , f'J, / e , f p , g p , g' p , g e , h p , and h' p , which do 



not depend on any kincmatical variables, plus the phases S® and S\, which are 
still functions of s n . 

It is worth noting at this point, that the phase shifts enter into the decay 
rate only through interference terms, which are proportional to sin</> or coscj). 
The azimuthal (j) asymmetry, which has an amplitude of only approximately 
10% of the the average (Fig. is principal source of the phase shift information 
as illustrated in Fig. ||. 

There exist three slightly different possibilities to extract the form factor 
and phase shift information from the data. Firstly, we can use the parameteri- 
zation given above in Eq. [| and perform individual fits in bins of s T respecting 
the s n dependence of the phase shifts. Secondly, there is the option to express 
the phase shifts by Eq. |J shown below. This allows to use the whole data 
sample in one single fit. Finally, we can employ in addition Eq. |5[ reducing the 
number of parameters by one. 

The phase shifts 6$ and S\ can be related to the s-wave scattering lengths 
<Jq and dp by means of the Roy equations HI) . A recent analysis of the phase 
shift data © used the parametrisation proposed by Schenk © : 



The Roy equations were solved numerically, expressing the parameters Ag, B],, 
Cj, D\, and s\ as a function of the scattering lengths <Zg and Oq. The possible 
values of the scattering lenghts were restricted to a band in the <2q — Qq plane, 



band, the universal curve can be expressed by a second degree polynominal, 
relating a[j and 



For the fits our data are divided over 6 bins in s T , 5 bins in s e , 10 bins in 
cos0 T , 6 bins in cos# e , and 16 bins in 0, 28'800 bins in total with on average 13 
events in each bin. In the \ 2 minimization procedure the number of measured 
events rij in each bin j is compared to the number of expected events in that 
bin rj given by 





a; 



2 




0.0849 + 0.232 a° - 0.0865 (a°) 2 



(5) 




(G) 



F MC , G MC , and H MC correspond to the original form factors used to gener- 
ate the Monte Carlo event, while F New , G New , and H New are the form factors 
derived from the parameters of the fit. q denotes the inital and not the recon- 
structed kinematical variables. Thus the parameters are applied on an event 
by event basis, and at the same time a possible bias caused by the initial choice 
of the matrix element is eliminated. 

In a first step the form factors in each bin in arc fitted independently. 
It is further assumed that the form factors do not depend on s e and that the 
form factor F contributes to s-waves only, i.e. f e — g e = f p — are fixed at 
zero. This leaves one parameter each to describe the form factors F, G, and H 
per bin plus the phase shift difference <5 = 6q — S{. Figures ^ and |^ summarize 
the results for F, G and the phase shifts, respectively. The dependence of F 
and G may be fitted with f s = 5.770±0.097, f' s = 0.95±0.55, f'J = -0.52±0.58 
(quadratic fit) or with f a = 5.835 ± 0.065 , f' s = 0.47 ± 0.15 (linear fit), and 
g p = 4.684 ± 0.092, g' p = 0.54 ± 0.20 (linear fit), respectively. 

If the individual phase shifts are expressed in terms of the scattering 
length <2q (Eq. |J), i.e. the universal curve is used to express a$ as a function 
of Oq, and the form factors are parametrized as in Eq. [|, the results listed in 
Table || are obtained. 

Table 2: Results of the fits to the K e ^ data using Eq. | and Eq. | ( X 2 /NdF = 
1.075j. For the last column the latter constraint is dropped. The other param- 
eters and the x 2 do not change significantly. 



f s 5.75 ±0.02 ±0.08 
f' s 1.06 ±0.10 ±0.40 
f'J -0.59 ±0.11 ±0.40 


g p 4.66 ±0.47 ±0.07 
g' p 0.67 ±0.10 ±0.04 
hp -2.95 ±0.19 ±0.20 


a" 0.228 ±0.012 ±0.003 (UB) - 0.036 ±0.009 


a% 0.203 ± 0.033 ±0.004 


-0.055 ± 0.023 ±0.003 



In conclusion it can be stated, that the E865 results constitute a major 
advance in the program of precision tests of hadronic interaction models at 
low energies. They also do nicely confirm the updated two-loop ChPT predic- 
tions 



10) 
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Figure 4: K e ^ control plots. From top left to bottom right: ttt: mass; cos 0„; 
ev mass; cos9 e ; missing v mass squared; <j>. Data (solid circles), Monte Carlo 
(histogram) including backgrund, and background separately (lower histogram) . 
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Figure 5: Form factors F and G as function of s. 




Figure 6: Left: E865 phase shift, results compared to the Geneva-Saclay 
data El) . The lower curve is the result of the fit to our data. Right: irir 
scattering length. The solid vertical line and the ellipse indicate our results for 
the scattering lengths for the two analyses requiring the universal curve (UB) 
as additional input or not. The solid squares indicate sucessive orders in chiral 
perturbation theory calculations, the latest of which 0) is marked with a cross. 
The solid circle indicates LQCD-result^^ . The dot-dashed lines indicate the 
center and lower error band of the previous experiment — ) . 



